QUANTUM CIRCUIT IDENTITIES 



CHRIS LOMONT 

Abstract. There are quantum circuit identities that simplify quantum cir- 
cuits, reducing the effort needed physically to implement them. This paper 
constructs all identities made from 3 or fewer operations taken from a common 
set of one qubit operations, and explains how they may be used to simplify 
the cost of constructing quantum circuit identities. 



1. Introduction 

Since the breakthrough algorithms of Grover |Hj and Shor their generaliza- 
tions to the hidden subgroup problem (HSP) by Kitaev and beyond llfij . 
there has been intense interest in building quantum computers. But building quan- 
tum circuits is a challenging task, and one important part is being able to reduce 
a circuit to simple building blocks that are easily implemented. Some blocks are 
more costly to implement in terms of time, effort, and cost, for different physi- 
cal realizations of quantum gates, so exchanging some operations for equivalent 
operations can be useful. Some different physical approaches are Raman-coupled 
low-energy states of trapped ions |S] , electron spins in quantum dots |S] , linear 
optics [J, nuclear spins in silicon |12|. and atomic cavity quantum electrodynamics 
[Jj. Furthermore, when building quantum simulators |llU15j . it is useful to simplify 
constructions to equivalent circuits using fewer or faster operations. A slightly old 
list of quantum simulators is at |f 8| . 

Often placing two larger circuit blocks in series allows the adjacent blocks to 
have some gates merged. This can be done by realizing when certain combinations 
of elementary operations can be replaced by equivalent operations which lower the 
cost metric. For a good introduction to quantum circuits see (SJ- 

Identities can be generated automatically with reasonable effort, but when au- 
tomatically reducing circuits, doing so at runtime is too slow for longer identities. 
It is therefore important to have on hand a large list of identities, so simplifications 
can be automated. The purpose of this paper is to start with a relatively flexible 
set of matrices, and construct all identities up to a given length from them. The 
value of this paper is a simple set of rules to simplify many quantum circuits au- 
tomatically and deterministically. This is done explicitly for one qubit identities 
made from products of 3 or less "elementary" gates. Since most simplifications in 
practice should come from merging a few gates in adjacent larger circuit blocks, 
choosing length 3 should encompass most needs. 
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This paper is laid out as follows: section El defines the matrices used in the iden- 
tities and their relevance, section explains the procedure to obtain the identities 
(filtering out trivial or redundant identities, and proving their validity), sectional 
contains a list of identities, and finally section [5] contains concluding remarks. 

2. Definitions 

Single qubit operations are in one-to-one correspondence with 2x2 unitary 
matrices, but certain of these matrices have special importance in the quantum 
computing literature O Chapter 4]. Most circuits in the literature are constructed 
from the following eleven unitary matrices: 

*-(; ■) J) -0 

R x (0) =cosf/-isin|X R y (6) = cos §7 - i sin §Y 

R z (9) = cos —7 — zsin 

X, y, and Z are the Pauli spin matrices, and R X1 R y , and R z are the correspond- 
ing exponential matrices, giving rotations on the Bloch sphere 2, Exercise 4.6]. H 
is the Hadamard matrix, used extensively in quantum algorithms. S is the phase 
gate, and T is called the tt/8 gate (even though n /4 appears in it). Trivially S = T 2 
and Z = S 2 ; they are each included since they appear in many quantum circuit 
constructions. P is a phase rotation by 9/2 to help match identities automatically. 

3. Procedure 

Infinite families of qubit operations are likely to be hard to construct physically, 
except in special cases, so we arbitrarily chose a finite subset of the matrices in 
section|21- small enough so that automatically generating identities is not too time 
consuming, yet large enough to yield a useful set of identities. We start with 35 
matrices A = {I,X,Y,Z, H, S, T, R x (6j ) , R y (9 3 ) , R z (9j ),P(9j)} where 9 3 - = 
for j = 1,2, ... ,7. Note that R x , R y , R z , and P all have period 47r and equal the 
identity when their argument is 0, thus the choices for 9j. Taking fewer choices for 
9j results in much fewer identities, but enlarging the number by a factor of 2 made 
the process unbearably slow. In order to shorten notation later, let Xj — R x (^-j), 
j = 1, 2, . . . , 7. Similarly define Yj, Zj, and Pj. 

The process to find the identities was programmed in both C++ and Mathemat- 
ica 4.2 to help weed out errors and avoid missing any identities. An identity is an 
equivalence Ihs — rhs, for example Z = HXH , where Ihs and rhs are strings of 
matrix names. The length of an identity is the number of matrices on the right hand 
side. Ihs always has length 1. To find all identities of length n, all n' A ' products of 
combinations of matrices from A are computed, and then the results are compared 
numerically with a small tolerance to each matrix in A. Numerical computations 
were a lot faster than doing accurate symbolic comparisons in Mathematica, and 
they were much easier to implement on the C++ side. This generated a list of 
identities of the form Ihs = rhs where Ihs is a single matrix, and rhs is a product 
of n matrices. Identities were passed through a filter (described in section l3~T|) to 
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remove redundant identities, trivial identities, and to clean them up by applying 
obvious simplifications. Finally, the identities of length up to 3 are listed in section 
modulo the filtering process. The final list of identities were verified symbolically 
in Mathematica 4.2 and numerically by the C++ program to ensure that filter- 
ing did not introduce errors. C++ source code and a Mathematica 4.2 notebook 
implementing the programs are at |14j , along with a file containing all (unfiltered) 
identities of length 4 and less, of which there are 400089 (see table GJ. 

3.1. Filtering identities. Each identity is run through the filtering process de- 
scribed below to remove redundant identities and those that can be shortened using 
previous identities. The filter performs the following steps on an identity: 

The Filter 

(1) Shrink: If the length d of rhs is 3 or larger, remove the identity if any 
identity with a length of 2, 3, . . . , d — 1 can be applied to shorten rhs. This 
way, knowing / = HH, we avoid identities like H = HHH which are 
deduced from shorter identities. 

(2) Negate: Apply the rules {Q 4 -> -I, Q 3 -> -Qj-a} for Q e {P, X, Y, Z} and 
j = 5,6,7. This replaces expressions like X 5 Y 6 — > (— -Xi)(— Y%) = X1Y2, 
making all subscripts less than 4, and shortening the list. 

(3) Clean: Merge the negative signs from the previous step, and move to the 
front of the string. This step is trivial mathematically, but must be per- 
formed symbolically in the computer. 

(4) Phase: Since the Pj commute with all matrices, merge them all to the front 
of rhs, except for a possible "— " in front. That is, replace expressions like 
P 2 = -PzSPxZx with I = -P 2 SZ 1 . 

(5) Normalize: Apply the 90 commuting rules {QQj — > QjQ, TZ — > ZT, SZ — > 
ZS, ST -> TS, SZj -> ZjS, TZj -> ZjT, VQ 4 -> Q^V] for Q e {X, Y, Z}, 
V G A, and the rules on identities {// — s- I, I to normalize the ex- 
pression, and to allow further simplifications. For example, this allows the 
simplifications ZZ3SZSZ2 — » ZZSZ3Z2 — * SZ$ — > —SZ\. 

(6) Collapse: Combine any obvious expressions using the 196 identities {QiQj - 
Q k } for Q e {P,X,Y, Z}, i,j £ {1,2, . . .7}, and k = (i + j) mod 8, where 
Qo = I- For example, X 2 X 3 becomes X 5 and Y 2 Y 3 Y 3 becomes /. 

(7) Merge: Remove trivial identities such as when the Ihs and rhs are symbol- 
ically identical, or Ihs = I. Remove duplicate identities at this step. 

(8) Rotate: To avoid a lot of identities, if the identity is now made up only 
of rotations Xj, Yj, and Zj, (and treating ±7 as rotations) it is discarded. 
Since these matrices act as rotations of the Bloch sphere they satisfy usual 
rotation identities, so this does not remove too many useful relations. For 
the online list they are retained for completeness. 

(9) Repeat: Repeat the above steps until the list of length n identities becomes 
stable. Repetition is necessary; for example step El may allow step El to 
simplify the identity further. 

(10) Grouping: After the list is stable, group similar identities to shorten the list, 
using symbols {A, B,C}, chosen cyclically from {X,Y,Z}. For example, 
the three identities I = XX, I = YY, and I = ZZ are replaced with the 
identity I = AA, and the three identities Z 2 = —XY, Y 2 = —ZX, and 
X 2 = —YZ are replaced with A 2 = —BC. 
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So roughly the filter applies known identities to try to shorten a given one, and 
returns standardized, simplified (using certain rules) identities. 

3.1.1. Filtering effects. Table^shows the effects of filtering on the number of iden- 
tities returned. The first row lists the counts with no filtering from section 13.11 
applied. In this case, both Mathematica 4.2 and C++ find 47 identities of length 1 
(almost all are trivial except X4 = Y4 = Z4 = P4 = —I). Going to length 2 added 
625 identities, giving the total 672 in the tabic, and length 3 added 15068 identities, 
totaling 15740. The C++ program was fast enough to find the length 4 identities, 
but the Mathematica program is too slow for this length and longer. This case 
needed to compute all |A| 4 = 35 4 = 1,500,625 combinations, and compare them. 
Length 5 and 6 are probably computable, but will take a lot of resources. 

The second row of the table is the number of identities returned when filtering is 
enabled, but not applying step |S] (rotation removal) and step 1101 (grouping similar 
patterns) from section 13. II The third line contains the counts when only step 1101 is 
skipped, and the last line has counts obtained by applying all filtering steps. 



Filter 


Length 1 


Length < 2 


Length < 3 


Length < 4 


No filtering 


47 


672 


15740 


400089 


Keep rots, No groups 


12 


66 


293 


1330 


Drop rots, No groups 


6 


54 


185 


982 


All filtering 


2 


36 


155 


931 



Table 1. Identity counts 



Looking at the table, and noting that even the 931 filtered length 4 identities 
are too tedious to put in a paper, I decided to put the more manageable 155 length 
3 and less identities. However the unfiltered lists can be found online [Tl). 

4. Identities 

The identities in table [5] arc all the identities resulting from 3 or fewer products 
of matrices from A in section [21 after filtering. They are sorted alphabetically on 
the right hand side for quick reference. When the symbols A, B, or C appear in the 
identities, that identity stands for the three identities where A, B, C are a cyclic 
permutation of X, Y, Z, as explained in step 1101 in section T3. II 

4.1. Hand simplification. Simplifying diagrams by hand is quite tedious, but can 
be assisted using these tables. Apply the procedure listed in section [3] to simplify 
the product until it simplifies no further, then look for identities in the list that 
apply. 

4.1.1. Example 1. We verify the identity Y = —XYX. Starting with XYX, we 
see no identities in the table starting with XY ... or YX . . . , so we look for a 
grouping pattern, and find A2 = —BC, which we apply as Z2 = —XY, giving 
XYX — > — Z2X. Expand using A2 = CB and collapse / = A A giving the complete 
transformation XYX -> -Z 2 X -> -YXX -> —Y, proving Y = -XYX. 
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J = 


AA 


Y = 


HY 3 Z 2 


z 3 = 


-P 2 ZxZ 


S = 


X 2 SY 2 


X = 


Y 3 HY 2 


/ = 


-AC 2 B 


Yx = 


HZ 


Z = 


P 2 Z 2 


Y = 


X 2 Y 3 H 


H = 


Y 3 HY 3 


A 1 = 


-BAxC 


Xi = 


HZ 1 H 


/ = 


P 2 Z 2 Z 


H = 


-X 2 Y 3 Y 


S = 


-Y 3 SXx 


Ax = 


-BA 3 B 


x 2 = 


HZ 2 H 


Zx = 


P 2 Z 3 Z 


H = 


~X 3 HZx 


H = 


-Y 3 X 


A 3 = 


BMC 


Y 3 = 


HZ 2 Y 


Zx = 


-P 3 S 


S = 


X 3 SY 3 


x 2 = 


Y 3 YH 


A 2 = 


-BC 


Y = 


-HZ 2 Y X 


Z 2 = 


-P 3 ZxS 


Y = 


X 3 YX 3 


H = 


Y 3 ZY 2 


A = 


BC 2 


x 3 = 


HZ 3 H 


z 3 = 


-P 3 Z 2 S 


Y = 


-X 3 ZXx 


Y = 


-Y 3 Z 2 H 


A = 


B 2 C 


z 3 = 


-PiZS 


I = 


P 3 Z 3 S 


Z = 


X 3 ZX 3 


H = 


Y 3 Z 2 Y 


A = 


B 3 AB 3 


s = 


P 1 Z 1 


s = 


P 3 Z 3 Z 


Yx = 


YHX 2 


Y 3 = 


-ZH 


A = 


B 3 CB X 


A 2 = 


-P 2 A 


A 4 = 


Pi 


Z 2 = 


—YHY 3 


Xx = 


ZXxY 


A 3 = 


-CA 3 B 


A 3 = 


-P 2 A X A 


I = 


PiA 4 


Y 3 = 


YHZ 2 


x 3 = 


-ZXxZ 


A 1 = 


-CA 3 C 


A = 


P 2 A 2 


x 3 = 


-SHS 


S = 


YSX 


Yx = 


-ZY 2 H 


A 2 = 


CB 


I = 


P 2 A 2 A 


z = 


SS 


x 3 = 


-YXxY 


I = 


-ZY 2 X 


A = 


~CB 2 


Ai = 


P 2 A 3 A 


x 2 = 


-sxs 


Y 3 = 


YX 2 H 


H = 


-ZY 3 


A = 


CxBC 3 


Y 3 = 


-P 2 HX 2 


H = 


SXxS 


H = 


YX 2 Y 3 


I = 


-ZY 3 H 


A = 


-C 2 B 


I = 


P 2 HX 2 Y\ 


X = 


SX 2 S 


I = 


-YX 2 Z 


H = 


~ZxHX 3 


A = 


C 3 AC 3 


x 2 = 


P 2 HY 3 


Y 2 = 


-SYS 


Zx = 


YZxX 


X = 


ZxXZx 


I = 


HH 


I = 


-P 2 HY 3 X 2 


Y = 


SY 2 S 


z 3 = 


-YZxY 


Y = 


-ZxXZ 3 


Y 3 = 


-HX 


Yi = 


P 2 HZ 2 


S = 


TT 


Zx = 


-YZ 3 Y 


Y = 


ZxYZx 


Zx = 


HXxH 


/ = 


-P 2 HZ 2 Y 3 


Yx = 


XH 


X = 


YxH 


H = 


-Z 2 HX 2 


z 2 = 


HX 2 H 


Y 2 = 


-P 2 Y 


H = 


XYx 


s = 


-YxSX 3 


Yx = 


Z 2 HY 


Y = 


-HX 2 Y 3 


Y 3 = 


-PiYxY 


Y 2 = 


-XY 3 H 


Y = 


-YxX 2 H 


Y = 


-Z 2 HY 3 


z 3 = 


HX 3 H 


H = 


P 2 YxZ 2 


z 3 = 


-XZxX 


z 2 = 


YxYH 


Yx = 


-Z 2 YH 


Yi = 


-HYX 2 


I = 


P 2 Y X Z 2 H 


Zx = 


-XZ 3 X 


H = 


-YxYX 2 


H = 


-Z 2 YxY 


Z = 


HY\ 


Y = 


P 2 Y 2 


z 3 = 


XZ 3 Y 


H = 


YxZ 


H = 


-Z 3 HXx 


Y% = 


—HY 2 H 


I = 


P 2 Y 2 Y 


H = 


-XxHZ 3 


H = 


Y 2 HY 2 


X = 


-Z 3 YZx 


Y x = 


-HY 2 X 


z 2 = 


P 2 Y 3 H 


S = 


XxSYx 


Z = 


Y 2 HY 3 


Y = 


Z 3 YZ 3 


X = 


-HY 3 


I = 


~P 2 Y 3 HZ 2 


Y = 


XxYXx 


S = 


-Y 2 SX 2 






Yx = 


—HY 3 H 


H = 


-P 2 Y 3 X 2 


Z = 


-XxYX 3 


H = 


Y 2 XY 3 






I = 


-HY 3 X 


I = 


~P 2 Y 3 X 2 H 


Z = 


XxZXx 


Z = 


—Y 3 H 






z 2 = 


HY 3 Y 


Yi = 


P 2 Y 3 Y 


Y 3 = 


X 2 HY 


Y = 


-Y 3 HX 2 






Y 2 = 


-HY 3 Z 


z 2 = 


-P 2 Z 


H = 


-X 2 HZ 2 


x 2 = 


—Y 3 HY 







Table 2. The 155 filtered identities up to length 3 



4.1.2. Example 2. Another way to approach this is to try to get all expressions 
into X,Y,Z,Xj,Yj,Zj, and then commute them to get X's together, etc., until 
simplified. Thus to simplify HXH , we have 

HXH g=XYl » XY 1 XXY 1 XY X Y X X Y 2 A=BC2 1 Z 



where the rules applied are above each arrow, resulting in Z = HXH . 
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5. Conclusions and future work 

There are several ways to speed up the search process, like adding the phase 
matrices only to the front, and not using the identity except for comparisons. But 
the speed improvements are minimal. 

Other directions are to extend this to understanding the 2 and 3 qubit operations. 
For example, the Toffoli gate can be implemented with 5 basic 2 qubit gates I 
believe it is unknown if this can be done with 4 gates, although it seems unlikely. 
A computer search should be able to shed light on this, and perhaps open up new 
understanding about minimal number of quantum gates needed for some other 
constructions. 

Finally, length 2, 3, and 4 filtered and unfiltered identities are online as weu 
as the C++ code and the Mathematica 4.2 code. 
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